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Abstract. We study an initial-boundary value problem (IBVP) for a coupled Cahn-Hilliard-Hele-Shaw 
system that models tumor growth. For large initial data with finite energy, we prove global (local resp.) 
existence, uniqueness, higher order spatial regularity and Gcvrey spatial regularity of strong solutions to 
the IBVP in 2D (3D resp.). Asymptotically in time, we show that the solution converges to a constant 
state exponentially fast as time tends to infinity under certain assumptions. 



1. Introduction 

Approximately 550,000 patients will die from cancer before the ball drops in New York City's Time's 
Square again. In the United Kingdom, one in four people will die of cancer, whilst one in three will at 
some point in their life be diagnosed to have cancer ( http: / / www.cancerresearchuk.com/) . In comparison 



to molecular biology, cell biology, and drug delivery research, mathematics has so far contributed relatively 
little to the area. A search in the PubMed bibliographic database (http: / / www.ncbi.nlm.nih.gov/PubMed/) 
shows that out of 1.5 million papers in the area of cancer research, approximately 5% are concerned with 
mathematical modeling. However, it is clear that mathematics could make a huge contribution to many 
areas of experimental cancer investigation since there is now a wealth of experimental data which re- 
quires systematic analysis. At the current stage of cancer research, most of the mathematical models 
are built and developed from three perspectives: discrete (microscopic), continuous (macroscopic) and 
hybrid (micro- macroscopic) . This paper focuses on a continuous mixture model. The model is designed 
to capture the dynamics of morphological changes in solid tumor growth. The thermodynamically con- 
sistent model is derived using a classical continuum mechanics approach, based on the principle of mass 
conservation and the second law of thermodynamics. The free energy is chosen by taking into account the 
effect of cell-cell and cell-matrix adhesion. In this approach, sharp tumor/host interfaces are replaced by 
narrow transition layers. These models are capable of describing the avascular, vascular and metastasis 
stages of solid tumor growth. In non-dimensional form, the simplest version of the equations reads (c.f. 

mm)- 

At := F'{(f>) - e'^Acj) ^ (j,^ - cj) - e^Ac/), 
V • V = S". 



(1.1) 



Here, (j) is the scalar volume fraction (order parameter) of tumor cells, /x is the chemical potential, S is 
a mass source term that accounts for cell proliferation, v denotes the advective vector velocity field, P 
is the scalar pressure, and e > 0,7 > are constants with the former describing the interface thickness 
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and the latter describing cell-cell adhesion. The function F{(j)), which appears in the definition of the 
chemical potential, i.e. in the second equation of p. II) takes the form F{(f>) = j{(t>^ — 1)^, which has a 
physically relevant, double-well structure, each of them representing the two phases of the mixture. We 
observe that, when 5 = 0, system reduces to a simplified version of an earlier model derived in 

[m [H] , that is used to describe gravity-driven, density- mismatched, two-phase flow in a Hele-Shaw cell. 
A similar set of equations was also used in [3D| as a model for spinodal decomposition of a binary fluid 
in a Hele-Shaw cell. In this paper, we are interested in the case in which 7 and e are strictly positive. 
Therefore, without loss of generality, we assume that e = 7 = 1 throughout the paper. 

We consider system in a bounded domain ft which is a rectangle in or a box in R'^. The 

system is supplemented by the following initial-boundary value conditions: 

r(^,A^)(x,O) = (0°,/)(x), xel]; 
(1.2) J v-n|aj2 = 0, 

[ V0 • n\on = V^i • n|ao =0, t>Q, 

where fl = [0,£] x [0,/] in and = [0,L] x [0,/] x [0,h] in R"^, and n is the unit outward normal to 
dQ. In view of the definition of fi and the boundary conditions p.2p we see that, in fact, the boundary 
condition for fi is equivalent to V(A0) • n|gn = 0. 

System (jl.ip is closely related to the Cahn-Hilliard-Navier-Stokes type system of equations studied in 
[aiSllllEllelllOlllSlEllllllMllMlESand many others. In the absence of the source term, system (fTTI) 
serves as a model for cell sorting owing to differential cell-to-cell adhesion [34] . 

The CHHS system (jl.ip . with 5 = 0, was studied numerically in |33j . where a convex-splitting scheme 
was proposed and was proven to be unconditionally stable and solvable. A practical and efficient multigrid 
method was used for solving the scheme at each time step. However, to the best of our knowledge, system 
(jl.ip has not been studied analytically, even when 5 = 0. In this paper, we will study qualitative behavior 
of solutions to (|l.ip - (|1.2p with 5 = 0, i.e., existence, uniqueness, regularity and long-time behavior of 
solutions to (|l.ip - (|1.2p with 5 = 0. Our first goal is to prove well-posedness of strong solutions, see 
Definition 2.1 below, to the IBVP. We show that for initial data belonging to H^, there exists a unique 
global (local resp.) strong solution to (|l.ip - (|1.2p of the 2D (3D resp.) system. The second goal of this 
paper is to improve the spatial regularity of the strong solution for the same initial data. Following ideas 
used for parabolic PDEs and in particular for the Navier-Stokes equations (c.f. we show that 
the spatial regularity of the solution increases by degree 2 within the lifespan of the solution. Thirdly, 
following [ini HZ] (see also [H]), we study the Gevrey regularity (spatial analyticity) of the solution in 
the domain fi. We show that for initial data belonging to H^, the solution belongs to the Gevrey space 
Gj (ri), see definition in Section 8, globally in 2D and locally in 3D. In other words, the solution is globally 
(locally resp.) spatially analytic in 2D (3D resp.). Since the elements in the Gevrey class have high-mode 
coefficients that decay exponentially in wave number to zero, as a result of the Gevrey regularity one can 
show that the standard Galerkin method, based on basis functions in ()2.ip - ()2.3p . converges exponentially 
fast, see e.g. [HI [121 EQJ [22l [27]. The last part of this paper is devoted to the study of the long-time 
dynamics of solutions to (|l.ip - (|1.2p . First, we show that cj) converges exponentially rapidly to its spatial 
average (f> over the domain O in the norm, as time goes to infinity, provided that cf) lies outside the 
spinodal region where F"{-) = 30^ — 1 > and the initial perturbation — (^||^2 is sufficiently small. 
The result holds true in both 2D and 3D. It should be pointed out that the idea for our proof is in the 
spirit of [3] , where the author studied the long-time behavior of strong solutions to an initial-boundary 
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value problem for a coupled Cahn-Hilliard-Navier-Stokes system. On the other hand, we show that, in 
the 2D case, without the smallness assumption on the initial perturbation, (f> still converges exponentially 
fast to (/) in the norm, as time goes to infinity, provided that I < L < ir. This condition indeed 
implies that the constant Co in Poincare's inequality (i.e., ||/ — /III2 < C'o||V/!|^2) on the domain 
n = [0, L] X [0, 1] satisfies Cq < 1 since, for the domain under consideration, it can be explicitly computed 
that Cn = (Ai)^^ = L^/tt^ using the eigenfunctions (|2.1[) given below, where Ai is the first eigenvalue of 
the Laplacian operator subject to the Neumann boundary condition. Regarding the condition I < L < n 
(or Cji < 1), we remark that, since the interface thickness e is assumed to be 1 in the beginning, the 
units arc still consistent. Indeed, in terms of e, the condition is I < L < en for the original model (jl.ip . 

Here, throughout this paper, || • \\lp, \\ ■ and || • denote the norm of the usual Lebesgue 

measurable function spaces {1 < p < 00), L°° and the usual Sobolev space , respectively. For 
p = 2, we denote the norm |j • ||^2 by || • || and || • \\w=-^ by || • \\h^, respectively. The function spaces under 
consideration are 

C{[Q,T]-H-m and L\[Q,T]- H'^n)) 

equipped with norms 

sup \\^{;t)\\Hr and ( / \\^{:T)\\l.dTY'\ 
0<t<T Jo 

where r, s are positive integers. 
Let 

V:={0e I \/9 ■ n\an = V(A0) • n\an = 0}, 

and 

H := Closure of V in the L" — norm, 
V := Closure of V in the — norm, 

W := Closure of V in the H'^ — norm. 

Throughout the paper, unless specified, C will denote a generic constant which is independent of the 
unknown functions ^ and v, but may depend on D,, initial data and the time T. The value of C may 
vary line by line according to the context. 

Our main results, for with 5 = 0, are stated in the following theorems. First, for the existence 
of strong solutions, we have 

Theorem 1.1 (2D Global existence). Let Q = [0, L] x [0,/] and let (jp ^ V he given. Then for every 
T > there exists a global strong solution {(j), fj.,v, P), see Definition 2.1 below, to the initial-boundary 
value problem m^-^T^ such that 4> € C{[0,T]]V) n L^{[0,T];W), /i € C{[0,T];H) n L^{[0,T];V), 
v£L^i[0,T];iH\n))^) andPeL^i[0,T];H^in)). 

Theorem 1.2 (3D Local existence). Let n = [0,L] x [0,/] x [0,h] and let e V be given. Then 
there exists a local strong solution {(f>, fijV, P) to the initial-boundary value problem il. l\) - n~^) such that 
(f> e C{[0,T*];V) n L\[0,T*];W), fi € C{[0,T*]; H) n L^[0,T*];V) , v e L^{[0,T*]; {H\n)f) and 
P G i^([0, T*]; i7^(fi)) for some finite T* > 0, which depends on the initial data and the domain ft. 

Concerning the uniqueness and well-posedness of the strong solutions obtained in the above theorems, 
we have 

Theorem 1.3 (2D & 3D Uniqueness and well-posedness). Let the conditions of Theorem 1.1 and The- 
orem 1.2 hold. Then the solution obtained in Theorem 1.1 (Theorem 1.2 resp.) is unique and depends 
continuously on the initial data. 
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For higher order spatial regularity of the solutions obtained in Theorem 1.1 and Theorem 1.2, we have 

Theorem 1.4 (2D & 3D Higher order spatial regularity). Let the conditions of Theorem 1.1 and Theorem 
1.2 he satisfied and let (j) he the strong solution to hl.l\ )- fT^) . Then 

m-MH^<^, v<e(o,r), 

f r'^+^U{.^r)\\]jUT<C{T)-. Vte{0,T), 
Jo a 

where C{T) is an increasing function ofT, a is an arhitrary positive constant, and T > is the lifespan 
of the solution, i.e., T > is arhitrary for the 2D case and 7" = T* for the 3D case. 

The next theorem gives the Gevrey regularity of the solution to (|l.ip - (|1.2p . 

Theorem 1.5 (2D & 3D Gevrey regularity). Let (ff^ he given. Then the solution to il.l] )- nT2\) is 

regular in the Gevrey sense, see definition in Section 8, glohally (locally resp.) in 2D (3D resp.), i.e., the 
solution is glohally (locally resp.) spatially analytic in 2D (3D resp.). 

For long time behavior of solutions to ()l.ip - ()1.2p . we have the following three theorems. 

Theorem 1.6 (2D Long time behavior for small initial perturbations). Let (jp he given and assume 
that (j) ~ 1^ J (fP{x)dx > \/3/3. Then there exists a constant S > 0, sufficiently small, such that if 
\\(fp — (^11^2 < S then the glohal solution to U.l\) - nT^) satisfies 

m-,t)-^l,<Ce-''', Vt>0, 

for some constant C > independent of time. 

Theorem 1.7 (3D Long time behavior for small initial perturbations). Let the condition of Theorem 1.6 
he satisfied. Then the strong solution to lll.l\) ~ !!T^) exists glohally in time and it satisfies 

m-,t)-^'jfi<Ce-''\ Vi >0, 

for some constant C > independent of time. 

Theorem 1.8 (2D Long time behavior for large initial perturbations). Let (jp £V be given and suppose 
that I < L < t:. Then the global solution to U. 1]) - ^^) satisfies 

Il0(-,t)-<^||?f2 <Ce-^*, Vt>0, 

for some constant C > independent of time. 

Remark 1.1. Our results show that the diffuse interface model is glohally (locally resp.) well-posed in 2D 
(3D resp.). Moreover, under certain conditions on the initial data, the solution collapses to a constant 
equilibrium state as time tends to infinity. This .suggests that the distinction between the tumor and the 
surrounding tissue in the microenvironment will blur as time proceeds. Theorems 1.6 and 1.7 indicate that 
constant solutions are locally asymptotically stable when they are outside the chemical spinodal. While, 
based on our analysis, it is unknown what happens if they are in the spinodal region where F"{(f)) < 0. 
This launches a new interesting problem for future pursue. In addition, we expect to extend the results 
to the case with S ^ by adopting the ideas in this paper. 
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We prove the above Theorem 1.1-Theorem 1.5 by combining the standard Galcrkin approximation and 
method of energy estimate. The energy estimate is dchcate mainly due to the coupling between the Cahn- 
Hilliard equation and the Helc-Shaw cell equations by advcction and order parameter 0. Great efforts have 
been made to simplify the proofs. The proofs of Theorem 1.1-Thcorcm 1.5 involve intensive applications 
of Sobolev and Gagliardo-Nirenberg type inequalities, see Lemma 2.2. The proofs of Theorem 1.6 and 
Theorem 1.7 are in the spirit of [1]. Since cj) is assumed to be outside the spinodal region (Jj) > VS/S), 
we note that F"{4>) = 3(/)^ - 1 > (c.f. (|l.ip 9') on a small neighborhood = - (5, + 5], where 5 > Q 
is sufficiently small. The idea is to solve a modified problem with F replaced by an auxiliary function 
whose second order derivative is non- negative on M and coincides with F"{-) on I^. Then, under 
the smallness assumption on the initial perturbations, it can be shown that the solution to the modified 
problem is indeed the solution to the original problem and converges exponentially to the constant state, 
as time goes to infinity. The last theorem is proved by careful exploration of the condition I < L < n and 
previously established a priori estimates. The proof involves exhaustive coupling of energy estimates. 

The plan of the rest of this paper is as follows. In Section 2, we give some basic facts that will be 
used in the proofs of Theorems 1.1-1.8. In Section 3, we present some lower order a priori estimates for 
the approximate solutions constructed through the Galerkin method. Then we prove Theorem 1.1 and 
Theorem 1.2 in Section 4 and Section 5 respectively. In Section 6, we show the uniqueness of the solution 
by defining a stored energy. Then we improve the spatial regularity of the solution and give a proof 
of Theorem 1.4 in Section 6. Section 7 is devoted to the proof of the Gevrey regularity of the solution 
following ideas from [151 HZl [29] . Finally, we show the long time asymptotic behavior of the solution in 
Section 9. 

2. Preliminaries and functional setting 

In this section, we give some basic functional settings and preliminaries for the problem. First we 
define the function space in which the approximate Galerkin solutions are established. For more details, 
we refer the readers to [27l [29] and references therein. For simplicity, we only give the 3D case, while the 
2D case follows similarly. We let 

H,n = span{(/)ij^fc|A,;j_fc < A,„}, 

where 



(2.1) 0,,,.(x) y 1^ cos (^x) cos (^y) cos 

for i,j, fc = 0, 1, 2, ... are the eigenfunctions of the Laplacian operator on the domain il = [0, L] x [0, /] x 
[0, h], subject to the boundary condition V0 • n = on dil with the corresponding eigenvalues 

/Z7r\2 /77r\2 /fc7r\2 

K,3,k = ( — I +(^) +(— ^ 



LJ \ I J \ h 

and < Ai < A2 < • • • < A„j < • • • denote the set of distinct eigenvalues Xij,kS ordered by their 
magnitude. 

Next, we define 

V := {v e {C°°(n))^\v ■ n\aa = and V • v = on 



and let 



H := Closure of V in the — norm, 
V := Closure of V in the — norm. 
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In addition, let 



(2.2) 



ik /ZTT \ /JTT \ /fcTT 

— sm — X COS —y cos — — 
L \L J \ I '^J \ h 



kir 



L 



y cos 



kn 



ir^ + p) '^'[tv ^°ht^J ^^Ht^ 



for i,j = 0, 1, 2, k = 1,2, ... and 



(2.3) 



7 /*7r \ /JTT \ /fcTT 

- sm — x COS —V cos — 
I \L J \ I "J \ h 

fin \ /JTT \ f kn \ 

{TV Htv ^°HT"J 





■ cos 



for i,j,k = 0,1,2,.... Then {vijys,\Vijys} form an orthogonal basis for H. Moreover, the finite- 
dimensional subspace is defined in the similar fashion as H„i, using the orthogonal basis {vj;j_fc's,Wi j-^^'s} 
(c.f. [ni[2ll[29]). 

Next, we let Po- : (L^(r2))'^ — > H be the Helmholtz-Leray orthogonal projection. By applying P^r to 
([OJa we have 

(2.4) v~P,{fiVcP) =0. 

The following lemma is given in |27| and references therein. 

Lemma 2.1. For given functions fi and (j), there exists a unique solution to the problem . I\1.1\ a 

subject to the boundary condition v • n|gf2 = 0. Moreover, there exists a constant C > such that 

||v||^.,, < C||AiV0||w-., s = 0,l, Vl<g<c». 

We also need the following Sobolev, Gagliardo-Nirenberg, Agmon and Poincarc type inequalities and 
interpolation inequalities which are standard and classic, c.f. IQ} [25 l I31j. 

Lemma 2.2. Let f = J^^ f dx. Then in the 3D case, there exists a constant C > such that 



(I) 
(11) 
(III) 

(IV) 
(V) 



.f\\h<C\\f\\l^, Vpe [1,6], 

f\\h<c{\\f\\'/^vfr^'+\\fr), 

f\\l^<C\\f\\%., 

f\\U<C\\f\\H4f\\H-, 

f\\h<c\\vf\\h, if 7 = 0. 



In the 2D case, there exists a constant C > such that 

(VI) !l/|li4<c(||/iu.||v/iu. + ||/||i0, 

(VII) ||/!li4 <q|/iu.||v/iu., if 7 = 0, 

(VIII) !l/||i=o <c||/iu.||/||ff2. 



ANALYSIS OF SOLID TUMOR GROWTH 



7 



In the 3D case, there exists a constant C > such that 

(IX) i/ilis <c(||/|U.||v/|U. + ||/||i.), 

(X) l/llis <q|/!U2i|v/!U2, if / = o, 

(XI) l/llioc, <c||/!U.j|/i|H3, 

(XII) I V/lli. <C||A/||i„ if Vf-n\an^O, 

(XIII) l/ll^. <C||A/||2,, if 7=0 and V/ • n|aj, = 0. 

Now we give the definition of strong solutions to the problem (|l.ip ~ (|1.2p . 

Definition 2.1. Let e V be given, and let T > 0. A strong solution of il. l]) ~ n~^) in the interval [0, T] 
is a pair of functions {(j), fi) such that € C{[0,T];V) n L'^{[0,T];W), fi e C{[0,T];H) n L^{[0,T]]V), 

^ e L^[0,T];H), and 

ot 

for almost every t £ [0,T], with fi given by (|1.1[) 9 and v given by ^2^. That is, for every lu £ H the 
above equation holds in the following sense 

(2.5) {^{t2),io) - {^{t,),io) + r (v- V(/.(T),c^)rfr- T (A//(t), w)dr, 

Jti Jti 

for every ti,t2 G [0,T], where ^ is given by ()l.ip 9 and v is given by \2.4^ . 

3. Lower order estimates in 2D and 3D 

In this section, we will give some general energy estimates which are valid in both 2D and 3D. The en- 
ergy estimates are performed for the approximate solutions to the original problem, which are constructed 
through the standard Galerkin procedure (c.f. [27 l I29j). 

We consider the Galerkin approximation system: 



(3.1) — "P- + Pto(v„ • V0™) = A^„, 

dt 

(3.2) ^rn = Prn{4>m) ~ 4'm - ^4>m, 

(3.3) V„i = P„iPff(^,„V(?!),„), 

with initial data 

(3.4) </)™(x,O) = 0°,(x) = P,„0O(x), 



where Vm G H™, (j)m,fJ-m & Hm, and Pm and Pm denote the orthogonal projections from H and H onto 
Hm and Hm, respectively. The finite-dimensional ODE system has a short time unique solution due to 
its locally Lipschitz nonlinearity. 

We now prove some lower order energy estimates which are valid in 2D and 3D. 

Step 1. Taking the inner product of p.ip with we have 

(3.5) ^(^^^ (i<^f„ - i0^)rfx+ i||V0™||2^ + W^flrnf = -(Pm(v™ • V</)™),/.„). 

By definition, it is straightforward to show that 

(3.6) - (f,„(v,„ • V0„O, /im) = -(v,„ • '^(t>m, IJ-m), 
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which, together with p.5p . imphes that 

(3.7) - ^€n)dx+ i|lV0™f ) + \\V^,„r = -(v™ . V^„,M. 
Next, taking the inner product of (j3.3p with v„i we have 

(3.8) ||v,„|P = (P,„P^(MmV0,„),V„,). 

Again, by definition we have 

(3.9) (P,„Pcr(AimV0„i), V„i) = {^J.mS7<j)m,V,n) = (v,„ • V(?!),„, ^m) , 

which yields 

(3.10) llVjiill^ = (v,ii • V(t)m,fJ-m)- 

Therefore, by adding p.7p and p.lOp we have 

(3.11) ^ - ^0™)rfx + i||V0™||2^ + ||Vm™||2 + ||v™||2 = 0, 
which is equivalent to 

(3.12) I (^lul - If + ^l|V0™ir) + ||V/i,„|p + ||v„|r = 0, 
and this in turns implies, by integrating w.r.t. time, that 

'^110™ - If + ^l|V0™f )(t) + 1^ (||VA.„,f + ||v™f )ds 



(3.13) 



^ll</'™(0) - If + ^l|v<^™(o)f ) < cin, 



where C(r2, ||0°||jLfi) denotes a constant which depends on ft and In the last step we used the 

fact that ||<^m(0)||L4 < C||0™(O)||//i < C||(/«°||_hi (see Lemma 2.2). As a resuU of ^J^, it follows that 
(|3.ip - (|3.4p exists globally in time, for positive time, because the nonlinearity is locally Lipschitz and the 
solution is bounded. 

Step 2. In this step, we derive some by-products of p.l3p . Integrating equation p.ip over 17 x [0,t] 
we have 



(3.14) / (/.,„(x,i)rfx= / (/.°„(x)dx= / (^"(x)dx. 

Jn Jn Jn 

Define 

(3.15) ^™ = <^=p^0°(x)dx. 
From p.l3p we easily see that 

(3.16) UJ\l^ = !|0™f +!|V0„J|2 < |O|i/2||0^„||2^+||v0„J|2 < C|!02^-l||2+j|V<^™||2 < C(17, 
which implies that 

(3.17) Umfm<C{n,U'>\\Hi). 
Then we have, by virtue of Lemma 2.2 (I) 



(3.18) 



</'L(x,t)dx 

n 



< \n\'^''Urn\\h < CU^nW'h < ^(17, ||0"|UO < C. 
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Thanks to the boundary condition V(/)m • n|ao — and p.l4[) 



(3.19) 



1 



/im(x,t)dx 



1 

m 



1 



, (x, t)(ix — (/). 



Then, from p.lSp we know that for any t > and m > 1. it holds that 

(3.20) |/2™| < C. 

Due to the boundary conditions and the trigonometric basis (j2.ip - (|2.3p . by using Poincare's incquahty 
one can easily show that 

(3.21) 110m - 4,,„\\h-2. < C||(-A)>,„||, ||0,„ - 0m||^2.+i < C||V(-A)^0„J, s > 1, 
and 

(3.22) f^rnU'^^ <C\\i-Ayfl„^\\, |1 Mm - A^m Ik^. + l < C| | V (- A)''m,„ ] | , S>1. 

Furthermore, if we write 0„i — <t>m — <t> and fim ~ Hm — Am, then p.ip - (|3.3p is equivalent to 



(3.23) 



dt 



+ Pmi^m ■ V0m) = A/im, 
Pm ((0m + 4>)^) - 0m - - A0m - p,„ 



= Pra [(0m + 0)' - (0m + 0)^] 
^ Vm = PmP<T(AmV0m)- 



A0„ 



Notice that in order to obtain the last equation we used the fact that fLm is x-independcnt and that 
PCT(V0m) = 0. From now on we will call 0m and flm by 0m and ^m, respectively, and use the above 
system for the next estimates, taking into considerations that 0m ~ 0, jlm ~ 0, and using the already 
established estimates p.l6p and p.20p . 

Step 3. Taking the inner product of (|3.23P i with 0m and using (|3.13p we have 

1 d 



(3.24) 



2dt 



||0m|P + 3|l(0m + 0)V0mf + ||A0m|P = ||V0m||' < ^(17,0°), 



where C(f2,0*') denotes a constant depending only on ft and initial data 0°, and we have used the fact 
that (-Pm(vm • V0m),0m) = 0. Integrating p.24p in time we have 

ft 



1 



(3.25) 



re[0,t] 



1, 



- sup |10m(T)f + / (3||(0m + 0)V0mf +||A0mf)ds<C(r!,0°)t+-||0m(O)|| 



< C(17,0")(t + 1). 



Especially, by ((3?2T|| and ([3?25)) we have 
(3.26) 



"m.||i/2 



%2ds < C. 



(3.27) 



Step 4. As consequences of p.l3p and p.26p we have the following. First, by (|3.23p 9 we have 

||VA0m||' = ||VPm((0m + 0)') " V0m - V//m||' 

< C(||V((0m + 0)')|P + ||V0m|P + llV/imlP) 

< C(||(0m + 0)'V0m||' + ||V0m|P + ||V//m||') 

< C(||0^V0m||' + ||V0m||' + IIVAimlP). 
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By Holder's inequality, Lemma 2.2 (I) and (|3.2ip we then have 

R.H.S. of < C{U,n\\U\\V<pm\f + ||V0™f + llV^mf ) 

< C(||V0™||2||A0,„||2||V0,„||2 + ||V0,„||2 + WVfMnf) 

fo 291] 

< C(||V(/),„f ||A0„,||2 + ||V0„,!p + ||Vm™!P) 

< C(||A0,„||2 + ||v^„ 



2^ 



where we used Lemma 2.2 (IV) and also applied (|3.13p for the estimate of ||V(/)„i||^. By plugging (|3.28p 
into p.27p we have 

(3.29) ||VA0,„|p < C(|jA0„f + llV^^f ), 
which, together with p.l3p and (|3.26p . yields 

(3.30) / \\<l>rn\\H.ds<C. 



We note that the constant C on the RHS of ()3.29|) is independent of time. 
Second, similar to p.27p - p.28p we have 



< C{\\V{cj)rn + ^ff + ||V0™f + II A0,„||2) 

<C(||(0™ + <^)V</)„||2 + ||A0,„i|2) 
<C^I|A0™f, 

where we used Lemma 2.2 (XII). It also holds that 

||VM™f < C{\\{4>.m + ^f^4>mf + ||V0™f + ||VA0„f) 

(3.32) <C(||A0„,||2 + j|VA</)™f) 

<q|VA,^,„||2. 

We remark that, according to the definition of strong solutions, with the estimates obtained above, it 
is still not enough to take limit, in m, and conclude the existence of solutions to the original problem due 
to the nonlinearity in the equations. Therefore, we have to seek higher order estimate of the approximate 
solutions. 

Step 5. Taking the inner product of p.23P i with A'^ipm and applying the Cauchy-Schwarz and 
Young inequalities we have 

(3.33) = - / (v,„- V(/)„O(A20„,)dx+ / A(((/.,„ + ^)3)(A2(/,„,)dx- / A0„(A20„)dx 

Jn Jn Jn 

<i|| A20„||2 + ^||(v„ . V0™)||2 + ^WAiicf,^ + 03))||2 + ^\\A4>.mf. 
After rearranging terms we update above estimate as 

(3.34) ^ll^'^-ll' + < 3|l(v™ • V0„O||2 + 3||A((,^,„ + 4>fW + 3i|A(/.,„f . 
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''ml 1^4] 



ll ■ 



For the second term on the RHS of p.34p . by a direct calculation we have 

l|A((0,„ + 0)3)||2 

(3.35) <C(||(0™ + 0)||io.||A</)„||2 ^ ||(^^ ^ ^^||2^||^^^^ 

<C{Urn\\U + 1)I|A(/.™||' + C{Urn\\l^ + l)||V0„ 

Using Lemma 2.2 (IV), ((Xn)) and we estimate the first term on the RHS of as 

C{U^\\U+l)\\Acf,„,\\' 

<C{Urn\\HA\(^rn\\H-+l)\\^4'rnf 

<C7(|lA0„Jp + l)||A<^„J|2 
<C(||A0™f + ||A0™||2). 
Similarly, by using Lemma 2.2 (II) we have 



(3.36) 



(3.37) 



<C{ 
<C{ 
<C{ 
<C{ 



ll'/'mllioc, + l)||V</.™||i4 

UrnWmUrnWm + l) (| | 1 1 1 1 A(/)„J ^ 

||A0„,|| + l)(i|A0„,j|-Vi|V0,,„||-^) 
||A0,„|| + 1)||A(/) 

m 1 1 

||A0„||4 + ||A0™||2). 



||V0„ 



Plugging (p:^ and (|X77)) into ([5351) we have 

(3.38) ||A((0„ + 0)3)112 < c(||A0„ 

So wc update (|3.34p as 



(3.39) 



dt 



||A0„|p + llA^^^lp < 3||(v„ • V0„)|p + C(|1A0„|14 + llA^^f). 



We remark that so far we have been focusing on the 3D case with all the estimates. Next, we shall 
deal with the first term on the RHS of p.39p . Due to differences between Sobolev embeddings in 2D and 
3D, we have to separate the arguments into the two and three space dimensions, which will be given in 
the next two sections. 

4. Global existence in 2D 

In this section we prove Theorem 1.1, i.e., global existence of strong solutions to (|l.ip - (|1.2p in 2D, 
utilizing of the a priori energy estimates established in previous section. 
First, by (jXT^ and ([5?^ we have 

||(v™ • V(/.™)||2 < i|v,„||2j|V0,„||i,c 

(4.1) < C||AimV0„,|n|V0,„||||V0,„||^f2 

< c||/i™||i4||v0™||i4||vA0„||. 

We estimate the RHS of (|iTT|) as follows. For the first term, using Lemma 2.2 (VII) and (|53T|) - (P3^ we 
have 

(4.2) 11^*^1114 < Ci|//,„||i|VMm|| < C||A0,„||||VA(/),„||. 

Since V(?i>m ■ ^\dn = 0, the second term is estimated, using Lemma 2.2 (VI) and (XII) and p.l3p . as 



(4.3) 



||V<^, 



mil £4 



<C||V0„J||A<^,„|| <C||A<^,„|| 
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Combining (|ITT|) . and wc have 

(4.4) • V0„OII' < q|A0„|n|VA0,„||2. 

Plugging (|4.4p into p.39p and applying Lemma 2.2 (V) to A(p,n wc have 



(4.5) 



i^ll A0„|p + i|| A20„||2 < C||VA(/.„j2||^^^j|2 ^ c(i|A<^„,f + l|A,^™f ) 



< C||VA(/.,„|n|^^^j|2 ^ c||vA0,„||2. 
Applying Gronwall's inequality to (|4.5I) and using (|3.30p we have 

(4.6) sup ||(/.™(t)||^2 + / \\(t),nfH,ds < C. 

Te[o,t] Jo 

As consequences of (|4.6p we have 

||A'm|lH2 < C||AAi,„|P 

< A((0„ + 0)3)||2 + II A0„||2 + II A2^„||2) 

(4.7) 

<C(||A(/.„J4 + ||A0„,||2 + ||A20„J2)^ hy^, 

<C(||A0„r + ||A20„||2), bydSSIl, 

which implies, by (|4.6p that 



(4.8) 



/ \\^^r,l\\H2ds < c. 
JO 



From (14. 4p and (|4.6p we see that 



(4.9) 







||(v™ • V</)™)||2ds < C / ||A0™||2||VA0™fds < C / ||VA0,„||2ds < C. 



ds < C. 



Combining (|48)) - (|49)) and (|3^ i we have 
(4.10) 





d(j)m 


Jo 


dt 



Moreover, we have 



(4.11) 



||v„i|||jids < / (|lAimV0m|p + ||V^„iV0 

rn\\ 1 1 ^ y^m 



< 



(llA*m||i=c||V(/) 

771 \ \ 



Iiv0,„lli4 



< C / ||/^™||2,,||A(/)„||2ds 



< c 



\fJ-m\\H2ds < C. 







d4)n 



dt 



')ds < C{t), yt>o. 



Collecting above estimates we have 

(4.12) sup {UUr)\\h-+\\t^rn{T)W')+ (|| V„ || + || 0„ || 
T6[0,t] Jo 

Now, we are ready to use these estimates and apply standard compactness theorems to pass to the 
limit, in to, in order to show the global existence result in the 2D case. Fix T > 0. In what follows, all 
arguments will be carried out within the interval [0,T] unless specified otherwise. First, we note that, by 
(|4.12p . the sequence {0™} is bounded in L°°([0, T]]H'^) and L'^{[0, T]]H^) respectively, then by the weak 
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compactness theorem, there exists a g L°°{[0, T]; iJ^)nL^([0, T]; H'^) and a subsequence {0rn'} of {4>m} 
such that {0m'} converges to in the weak* and weak topology of L°°{[0,T]; H^) and L^([0, T]; 7?^), 
respectively. Then we have from (|4.12l) that 



(4.13) 

and 

(4.14) 



(s)||_H"2 < liminf ||0„i(.s)||//2 < C, for a.e. s e [0, i], 



ds < liminf 



\4'm{s)\\mds < C. 



In addition, since H is compactly embedded in H^^ and since the sequence {-^} is bounded in 
i^([0, T]; L^), by Aubin's Compactness Theorem we know that {(j)m'} converges strongly to (j) in i^([0, T]; i?"^). 
From now on, we denote all the subsequences of {4>m} by itself by extracting and relabeling. Thus, {4>m} 
also converges strongly to (j) in H'^ almost everywhere in [0, T]. In particular, 



(4.15) 



1 0m (t) 11/^3 — !► ||0(<)||^3 pointwise everywhere on i5 C [0,T]; 
\[0,T]\E\^0. 



By virtue of the estimate (|4.12p we know that there exists a function v G L^([0,T]; V) such that {vm} 
converges weakly to v in L^([0,T]; V). Therefore, we have 



(4.16) 

In summary, we have 



/ ||v(s)||^ds < liminf / 
Jo Jo 



A.s)\\U.s < C. 



(4.17) 



dt 
v„ 



weakly in L^{[0,T]; H'^), 
strongly in L^{[0,T]; H^), 
weak* in L°° {[0,T]; H^), 

^ weakly in ^^([o^ y]; ^2^^ 
at 

V weakly in L^{[0,T];Y), 



respectively. 

Since the sequence {0m} is bounded in I/°°([0, T]; H"^) and {^^f-} is bounded in £^([0, T]; L^), we know 
from the Arzela-Ascoli Theorem that {4>m} converges to in C([0,r];L^), and hence in Cvi/([0, T]; L^), 
Cvy([0,T];i?i) and Cw{[0,T]; H"^). In particular, we have 

(4.18) (0m(t2),C^) - (0m(tl),w) (0(^2), C^) " (0(il),w), 

for every ti,<2 € [0,7"] and every uj Cz H. Moreover, from the first weak convergence in (|4.17p we know 
that for every ti,t2 G [0, T] and every uj G H, 



t2 



(4.19) 

Next, we will show that 
(4.20) 



(-A^0m - A0m, W)ds 



t2 



(APm(0^),w)ds- 



(-A^0- A0,a;)ds. 



(A(03),u;)rfs-^O, 
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as m OD for every ti,t2 G [0, T] and every uj G H. First, we notice that, by Lemma 2.2 (III) and (14.17^ ^ 

ll0?„-0'll = ll(0™-'/')(0™+'/'m'/' + 0')ll, 

<2(||0„||2^+||0||2^)||0™-0|| 

(4 21) 

< C\\(t)m - 011- 

Second, by a simple calculation and (|4.12p one can show that 

< C(||(/)„,||ioo||A2,/,„J + ||0,„|U^||V0,„|U4||VA0,„|U4 

+ ||V0„ji6||A0„,||iC + ||0„,||ioc||A0„||Loo||A0„||) 

(4.22) < (:;(||(/)„J^2||A20,„|| + ||0,„||^2||V0,„||hi||VA0„Jhi 

+ l|V(/)m||^i||A0„,||^i + ||0„||h2||A0„||^2||A0„||) 
<C||A20,„||. 

And the same is true for cj), i.e., ||A^(0'^)|| < C||A^(/)||. Then we have 



||AP„(0^) - A(03)j|2 = (p„(^^) _ q^^,A'{P^{(bl,) ~ 0^)) 
(4-23) < \\P^{cl>l) - 0^|U2||A2(P„(0^) - 03)1,^^ 

< C(||0,„ - + ||(P,„ - I)q}^\\L2){\\A^cj),„\\L2 + |1A20|1^2 
where we have used (|4.2ip and (|4.22p . Integrating (|4.23p w.r.t. time we have 



t2 



tl 



||AP„(0^) - A(03)||2d5 < C(|10™ - 011^2 + ||(P„ - I)cb^\\L^) / (||A20„||i2 + ||A20||i2)rf, 



tl 



< 



C(|i(/',„-0||L2 + ||(P,„-/)03|,^^)^ 



which, together with the strong convergence of (f>,n in C([0, P]; L^) and the convergence || {P„i — I)<f>^ 
0, implies that 

(4.24) 



|L2 



lim 

m— ^00 



t2 



\AP^{(j>l) ~ A{cj)^)f-ds ^ 0, \fh,h e [o,r]. 



tl 



By Holder's inequality we have 
r-t2 



(4.25) 



{APraiK),^)ds~ / (A(03),C^)d5 



tl 



t2 



< 



|AP™(</>3j-A(03)||2d,)'ll^l|rt. 



Hence, by passing to the limit in m and using (j4.24p we have 



(4.26) 



lim 

m— >oo 



(AP„((/)j;,),L.)ds- / (A(03),w)ds 



tl 



t2 



0, Vii,t2 e [0,r]. 



This gives (|4.20p . As a consequence of (|4.19p and (|4.20p we then have 



(4.27) 

as m — >■ 00, where /i 



t2 



{Afim,io)ds— / {Aii,u})ds — > 0, 



A(j). 



Upon taking the inner product of p.23P i with w e iJ and integrating over the interval [ii, ^2] C [0, T] 
we obtain 



(4.28) 



tl 



i(t2),a;) - (0m(ti),a;) + / (P™(v„ • V^™), w)(is 



t2 



(A^m, w)(is. 
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In view of (|4.18p and (|4.27p we see that it remains to show 



(4.29) / (P™(v™ • V0™),L^)ds- / ((v-V</)),w)ds-^0, 

as m -> oo for every ti,t2 G [0,r] and every w e iJ, in order to prove (|2.ip . First, we notice that for 
every uj Cz H, 



(4.30) 







\\{W(p)uj\\'^ds < I ||V0||ioo||a;|pds 
/o 



^ads < C. 



< C\\Luf 

This imphes that (V(/))aj G L^([0,T]; (i^(f2))^). Moreover, by Holder's inequahty we have 



{P,n{'Vvi-'^(f>m),^)ds- / ((v-V0),w)ds 



(^'m(Vm • V(/)m) - Pm(Vm • V(^), w)ds + / (Pm(v,„ • V0) - Pm(v • V(/)), w)(is 



< 



+ / (P,„(v-V0)-(v-V0),L^)ds 



||v™-V(0™-</>)|||ltj|lds+ / ||(v„- v) • V</>||||tj||ds 



+ / ||P™(v-V(^)-(v-V(/))||||w||ds 



(4.31) 



<il^!l 



|v™|Ui|lV 

t2 



-</>)||Hic?s + ||a;|| / ||v,„ - v|l!|V0|lL-ds 

Pra{^-Vcj))-iv-Vcj))\\ds 



ids 



ds 



Wds 



\'\j3ds 



||i',n(v • V0) - (V • Vcj>)fdsj \h - i2| = . 

Since Vm is bounded in i^([0, T]; V) and the sequence {0m} converges strongly to in L^{[Q, T]; iJ"^), we 
know that the first term on the right hand side of (|4.3ip vanishes as m H> oo. For the second term, since 
Vm converges strongly to v in L^([0,r];H) C i^([0,r]; (L^(ri))^), we know that it also goes to zero as 
m — >■ oo. Furthermore, the last term tends to zero as to — >■ oo by virtue of the pointwise convergence of 
the projection operator. This proves (|4.29p . Thus by letting m oo in (|4.28p and using (|4.18p . (|4.27p 
and (j429l) we get 



(4.32) 



(0(<2),w)-(</'(tl),^)+ / ((V- V0),C^)ds 



(A/i, uj)ds. 



Moreover, by using the strong convergence of {0m} in L?' {[Q ^T]\ H^) and the estimate (|4.12p one can 
show that Vm ~ PmP<T(Mm^0m) also converges weakly to Pcr(/iV0) in L^([0,r];H). This gives that 
v = P<,(AiV0). 
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In order to complete the proof, it remains to show that E C([0, T]; iJ^). For this purpose, we first 
re- visit the estimate (|4.5p . which reads 



(4.33) 1 ^||A0,„||2 + i||A2^„f < C\\VAcj)^f\\Acj)^f + C||VA0„,||2. 

We observe that, by virtue of the boundary condition V0m • n\dn = 

|1VA(/.„||2 ^ _(A0^, A^c/.^) < |1A0„|1||A2(/,„||, 
which imphcs, by Young's inequahty, that 

(4.34) 

< -\\A'4>mf + C{\\A^„,\f^ + \\A<f>„M^). 

After rearranging terms we have 

(4.35) 2^||A0„f + ||A20„||2 < C{\\Ac^,n\f + ||A0„f ) < C, 

where, in the last step, we used (|4.12p . Integrating ()4.35p over [s,t] C [0,T] for s < t we have 

(4.36) + \\A^4>rr.{T)fdT < || A0,„(s)f + C(t - s). 

Due to (|4.15p and the inequality, 

liminf(ari) + liminf(6„) < limsup(an + 6„) 
for any sequences a„, 6„ > 0, we have 

(4.37) \\Am\' + f \\A^<t>{r)fdT < || A0(s)||2 + C{t - s), 

J S 

for every s,t G E, s < E is given in ()4.15p . After rearranging terms in (|4.37p we have 

(4.38) \\\A<j>{t)f-\\A<f>is)f\< f \\A^{T)fdT + C{t~s), 



which can be extended to all s and t in [0, T] by manipulating sequences in E that converge to elements in 
[0,T]\£:. Using (|4:38| and the fact that cj) € Ciy ([0, T]; i?^) wc get € C([0, T]; i?^). Finally, we notice 
that T > is arbitrary. This gives the globally existence in 2D. This completes the proof of Theorem 
1.1. 

5. Local existence in 3D 

In this section we prove Theorem 1.2. Now we turn to p.33p in 3D. Instead of applying the Cauchy- 
Schwarz inequality to all the terms on the RHS of p.33p . we keep the first term and apply the Cauchy- 
Schwarz and Young's inequalities to the other two terms to get 

i||lA0„f + l|A2,/,,„f 

(5.1) = - / (v,„ • V(/.™)(A20„)dx+ / A((0„ +^)^)(A2</.„,)dx+ / A ((/.„, )(A2(/.,„)dx 
Jn Jn Jn 

<hA''(bmf~ f (v„.V0™)(A20„)dx + 2||A((0„ + 0)3)f+2||A0™|p. 



■4' 



n 
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By applying Lemma 2.2 (IX), we estimate the second term on the RHS of (|5.ip as: 
(vm • V(/),„)(A^(/),„)dx 



<il(v™- V0™)||||A2,^„|| 

<l|Vm!lL3||V0,„|U6||A20„J 

<C(i|v,„||i/2||Vv,„||i/2 + i|v,„||)||V(/.,„||Hi||A2<^,„|| 

<C||/i™V(/.™||i/2(||V/^„ . V(^„.||i/' + ||Ai™V20„.||^/2 ^ ||^„V0™i|i/2)||A(^™||||A2</,™||, 

where in the last inequality we used p.23p ?{ and Lemma 2.2 (XII). 

For the first term on the RHS of (fO]) . by Lemma 2.2 (X) and (XII) and (|3?3T|) we have 

||/i™V0,„||l/2 <||Mm|li('||V0,„||^// 

<ciiA^„.ir/iivA.™ir/iiv0™ii^/f 

(5.3) 

<q|A0„,||i/^||AM™||i/^||A,^„||i/2 

<q|AM™||i/^l|A0„||3/^ 

In order to estimate ||A^„i||^/'', by (|4.7p we have 

(5.4) ||Aa.„J1/4 < C(||A0„,!|i/2 + l|A20,„||i/4). 
Plugging (|5.4p into (|5.3p we have 

(5.5) ||a^,„V</)„J1/2 < c(||A0,„l|i/2 + l|A20,„||i/4)j|A^,„||3/4^ 

Next, to estimate the three terms in the bracket on the RHS of (|5.2p . we use similar arguments and 
Lemma 2.2 (XII) and (XIII) to obtain 



(5.6) 



■rn II ^3 II V ^Jm\\ 1^6 
< |lV^,„||^f ||V(/),„||^^i^ 

<q|AM™||i/'||A0,„||i/2 

<C(||A0,„|| + ||A2,^,„||i/2)||A,^,„||i/2, 



liM™V2</>™||^/' < ||^™||^/i||A0„||i/2 



<qiA*,„ii^4'||A0„,||i/2 
<q|AM„,l|i/2||A(/.„ji/2 

<C(||A(/.„J+||A20„jl/2)||A0„||l/2^ 



(5.7) 



and 



(5.8) 



< \\fJ-rn\\H?\\y(t)'m\\Hl 

<C||A^,„||i/2||A0™||i/2 
<C(||A0™|| + ||A20™||i/2)||A0„||i/2. 

Therefore, we have 

(5.9) ||V/i™ • V</)„Ji/2 ^ ||^,„v2(/.,„||i/2 + ||^„,V(/.,„||i/2 < (||A(/),„1| + ||A2,^,„||1/2)||a,^,„||1/2. 



(5.10) 
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Combining (|5.5p and (|5.9p wc get 

||At,„V0™||l/2(||VM,n • V0,„||l/2 + ||MmV20,„||l/2 ^ || V0™ 1 1 ^Z^) 
<C(||A,/),„||3/2 + ||A2</),„||3/4)j|A0,„f /4, 

which, together with (|5.2p . implies that 

(5.11) - / (v,„ • V0,„)(A2(/,„,)dx < C||A0„j9/4||A20,„f/4 + Ci|A,^,„||i5/4||A20„||. 

The estimate for the second term on the RHS of (|5.ip foUows from (|3.38l) : 

(5.12) ||A((0,„ + 0)3)f <c(||A0„,f + ||A0„f). 

Plugging the estimates (|5.1ip and (|5.12p into (|5.ip and applying Young's inequality we get 



2di 



||A0„ 



lA^ 



(5.13) 

which implies that 
(5.14) 



<i||A2,/,„f + C||A0„|r/4||A20,„||V4+ 

C||A0,„||i5/4||A20,„|| + C(||A(/>„J4 + ||A0,„||2) 
<i||A20„j2 ^ C(||A(/.,„||i« + ||A0,„||2), 



-||A(/)™||2 + ||A20,,„||2 < C(||A0„,||i« + ||A0„,|| 



We remark that the constant C on the RHS of (|5.14p is independent of time. 
Let 



(5.15) 

Then from (|5.14p we have 
(5.16) 

which implies immediately that 
(5.17) y,n{t) < 

Thus for any m we have 
(5.18) 



y,n{t) = !|A<^,„(<)i|2 + l. 

2/m(0) 



(l-8y™(0)8Ct)i/8' 



V < t < 



82/™(0)8C 



(1/2)1/8 

since y„(O)<y(O) = |lA0"||2 + l. 

As consequences of above estimates we have 



16(||A(/)0||2 + 1)8C 



(5.19) sup {UmiT)\\ (II VmllfJ-i + ||0m||H4 + ||/im||ff2 + 



re[0,T* 



T* 



')ds < C. 



These estimates, combined with the compactness arguments given in the preceding section, imply the 
local existence result in 3D. This completes the proof of Theorem 1.2. 
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6. Uniqueness and continuous dependence in initial data 



In this section, we prove the uniqueness and continuous dependence on initial data of the solutions in 
both 2D and 3D. Suppose that there are two solutions, namely, (0i, /xi, Vi, Pi) and (02, M27 V2, P2) on a 
joint interval of existence [0, T]. We use / to denote the difference of two functions. Then the differences 
satisfy 



A = 0? - 02 



t + V • V(/)i + V2 • V0 = A/i, 
- A0, 
V = -VP + flV(l)i + 
V0 • n|ao = V/i • n\on = v • n\sn = 0, 
0(x,O) = 0°(x)-0O(x). 



Taking the L'^ inner products of 
we have 



with yu and (|6.ip q with v, respectively, and adding the results 



(6.2) 



,-,112 



/X2V0 • vdx 



= / /i2V0 • vdx 



/iv2 • S/(j)dx. 



6v2 • V/idx 



Define 



G(0) 



((02+e)'-0i)de. 



(6.3) 

Then it is easy to see that 

(6.4) G(0)>O for^eM, G(0) = iff = 0, and G'(0) = (0? - 0^). 

Since, by the existence theorems, ^ € -^^([0, P]; H), then € G([0, P]; if) and is an absolutely continuous 
function. Moreover, since G(0) is a polynomial, it is differentiable, and G(0) is absolutely continuous in 
time with value in H and ^^^^^ ^ 

ED) 

(6.5) 



= G'(0)^ e L2([0,r];iJ). Furthermore, one can also show that (c.f. 



G(0)dx= / G'(0)-^dx. 



Then plugging ()6.5p into (j6.2p and using ()6.4p we have 



G(0)dx 



VAir < ll/^2V<; 



IVAII 



I0V2I 



We estimate the two nonlinear terms on the RHS of (|6.6p as: 



(6.7) 



IM2V0II' 



;!'V2||"' < ||/-f2| 



Iie||v2||i3 



<C'||M2ipH2||V0||2 + G(||V0||^ 



|V2||hi, 



where we have used Poincarc's inequality since 0i and 02 have the same average over fi. Then we get 
-f-l 



(6-8) 4(il|V0ir-^ll0f + /^G(0)dxj+J||^|p + ^||VA|P<G(||M2||?^. + ||v2||^O(l|V0f + 1^^ 
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Next, taking the inner product of (|6.ip i with we have 
2di' 



(v • V0)0idx + / //A0dx 



(6.9) 



(v-V0)0idx-||A0|p + ||V0|| 



bl +0102 +02)M0rfx 



< ||^||||0i!|l~!|V0!| - ||A0||2 + + n^J + ^^^^ + 0^||l^||0|||1A0|| 



<i||^|p-i||A0|p + C(||V0|p + ||0|p), 



which gives 
(6.10) 

Combining 



|||0ir + i!|A0!p<i||^iP + C(||V0|p + ||0ir). 



(6.11) 



d /I 



dt\2 



iiv0ir 



and (|6.10p and since G(0) > we have 

+ lur + ^ G(0)dx) + ^ii^ip + Jii v/i|p < c(t) (^11 V0IP + m'+l^ G(0)dx) . 



where 

(6.12) C'(t) =C(i|/i2l|^2 + ||v2||ffi +1) satisfying / C{T)ds < C, Vte[0,r]. 

Jo 

Applying Gronwall's inequahty to (|6.1ip we get immediately that for any t £ [0, T], 

G{4>{t))dx 

'1, 



(6.13) 



^iiv0(t)ir + iii0(t)ip 



< cxp 



ni. 

a0||2 



1„ ~n 



c{T)ds\{-\\^n' + 7^wr+ G{4>'')d^ 



2" " " 2' 

Therefore, the solutions depend continuously on the initial data in the sense of (|6.13p . In particular, if 
0" = we obtain the uniqueness of the solutions. This completes the proof of Theorem 1.3. 

7. Higher order spatial regularity 

In this section, we show that the unique solution obtained in previous sections enjoys higher order 
spatial regularity. The idea is adopted from the proof for the Navier-Stokes equations in I5T|. 
Taking the inner product of (|3.ip i with A'^0m we have 



i|l|VA0„f 



(7.1) 



<-||VA2 
-2" 



VA20„||2 

A30„,(P,„(v,„ • V0,„) - A(0„ + 0)3 + A0„)dx 

V(P™(v™-V0„)-A(0, 
3.„„. , „2 



= - / VA 
Jn 



|2 + ^(j|VA0„ 



f + A<j>,n)dx 
3)f + |lV(v„.V0™)||2) 



||VA((0„- 

We break the estimate of the RHS of ()7.ip into several steps. 

Step 1. For the third term on the RHS of (|7.ip . after direct calculations we have 

|iVA((0„, + 0)3)f 



(7.2) 



< 



C(|lV0,„||ia + ||(0,„ + 0)|lio.||V0„||2,||V20„j2^ + |l(0,„ + 0)||i^||VA0„ 
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Using Sobolcv embedding, Lemma 2.2 (XII) and (|3.21l) we estimate the first term on the RHS of (|7.2p as 
(7.3) \\^cf>m\\le < C\\W^,n\\%i < C\\A^„^\\^ < C7|l A0,„ f || V || ^ . 

Similarly, for the second and third terms, using the triangle inequality, Sobolev embeddings. Lemma 2.2 
(XII) and (XIII) we have 

c(!l(0™ + 0)||i=o||v</)™||i4||v2,/,„l|2, + ||(0„ + 0)||4^llvA0™||2) 

<C{Urn\\h\\^'l^rn\\Hi\\^^^r,r\\h + 1 1 V</)™ 1 1 || V^^^ 1 1 + 1 1 ,^,„ 1 1 1 1 V A0,„ || ^ + ||VA0„,||2) 

(7.4) 

<C(||A0™||2||A</)„|p||VA0„,||2 + ||A0™||2||VA0™||2 + i|A,^™f ||VA0„,||2 + ||VA0,„|p) 
<C(||A<^,„||4 + l)||VA<^,„||2. 
Plugging ((73)) and (fTil into ([7^ we have 

(7.5) |!VA((0™ + < C{\\A4>.m\\^ + l)||VA<^™f . 

Using previous estimates (|4.12p and (|5.19p we obtain for both 2D and 3D, within the interval of existence, 

(7.6) ||VA(0fJ|l2<q|VA</>™|p. 
Step 2. To estimate ||V(vm ■ V(/)m)|P, first, we have 

(7.7) ||V(V,„ • V0„O||2 < ||v,„ • V^(bmf + ||V(V,„) • V(bmf. 

The first term on the RHS of (|7.7I) is estimated as 

(7.8) <C\\^^rn\\h\\^^^ra\\m 

For the second term, using Lemma 2.2 (XII) and (XIII) and p.2ip and (|3.22p we have 



||V(v™) • V0™||2 < ||V(v„)|n|V0™f 



L° 

< C||V(Ai,„V</.™)|l2||V(/)„,||^2 

< C{\\VfIrn ■ ^(b^mf + 1 1 Mm 0™ | P) 1 1 V A0™ 1 1 ^ 

(7.9) < C{\\VfIrn\\l4^(^rn\\l4 + 1 1 /i,„ 1 1 lo. 1 1 V^^^ 1 1 2) 1 1 V A0,„ || ^ 

< C{\\VfIrn\\m\\y(^rn\\m + 1 1 Mm 1 1 ^2 1 1 A</)„ 1 1 2) 1 1 V A0,„ | ^ 

< C||AMm|P||A0,„||2||VA(/),„f 

< C||AMm|P||VA</)„||2, 

where we have used the uniform estimate of ||A0„i|p. Plugging (|7.8p ~ (|7.9p into (|7.7p we then have 

(7.10) |iV(v,„ • V(brn)f < C(||v,„||^i + 1|Am,„!|')||VA0,„||2. 
Therefore, combining (|7.1I) . (j7.6p and (|7.10p we have 

(7.11) ^I|VA0„J2 ^ ||VA20„J2 < c(l + ||v,„||2^, + ||AMm||')||VA0„||2. 

Step 3. Taking the inner product of p.23P i with A^(/)„j and using the Cauchy-Schwarz inequality 
we have 

(7.12) il||A20,„|p + l\\A^„-,\\' < (ll A2</)„J2 + 11 A2((</,„ + 0)3)||2 + II A(v,„ . V0„,)||- 
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By a direct calculation, Sobolcv embcddings and Lemma 2.2 (XII) and (XIII) we have 

l!A2((</.™ + ^)3)|l^ 
<C(||V0™||i4||A0„J|i^ + ||(,^,„ + 0)||i^||A0,„||i^||A0„,||2+ 

|l(0„, + 0)||i^||V0„||i3||VA0™||i6 + ||(0„ + 0)111^ ||A2,/,„||2) 
^^■^^^ < C(||A0„,||4||A20,„||2 + ||A0,„||2||A2(/,„J2||A0„j2 ^ ||A20„j2||^^^^j|2^ 

|iA0„j2||A^^^^||2||^2^^^j2 ^ ||A0„J2||^2^^^J2 ^ 1 1 A0,„ 1 1 ^ 1 1 A^^^ 1 1 ^ + ||A20,„||2) 

<C(|1A0„||2 + 1)'||A20,,„||2, 
which implies, within the interval of existence, that 
(7.14) ||A2((0„ + 0)3)||2<c|jA20„f, 

where we used (|4.12p and (|5.19p for the estimate of ||A(/)„i|p. So we update (|7.12p as 
(7-15) il||A20„||2 + ^il A30„J2 < c||A20„j2 ^ 3|j^(v,„ . V0,„)f . 

Step 4. Next, we estimate ||A(v„i • S/(f>m)\\'^. First, we have 

(7.16) ||A(v„ • V0„O|p < 4(||v,„ • VA(/.,„||2 + \\V{w„,)A(bmf + || A(v™)V</)™|P). 
The first term on the RHS is estimated, using Lemma 2.2 (XII) and (XIII), as 

(7.17) ||v,„- VA(/.„,f < ||v,„||23||VA<^,„||2« < C||v„J^i||VA0,„|l^i < C\\^„,\\]j,\\A^(j)„,\\\ 

Similarly, we have for the second term 

|lV(v„)A0„|l2 < ||V(v„,)in|A(/)„,|lio„ 

< C{\\Vfl^ ■ V(t>,n\\^ + ||MmA0m|P)||A0m||lf2 

(7.18) < C(||V/i„,||i4||V^,„||i4 + !|Ai,n||i=.||A(/),„||2)||A2^,„||2 

< C(||Am™|P||A0,„||2 + j|AAi„,|n|A0„,||2)||A20„f 

<C\\A^i,n\\^\\A^cf,„,\\\ 

where we used the estimate of ||A(/)m||2. 

The estimate of the third term on the RHS of (|7.16p is more delicate than previous two estimates. 
First, we have 

||A(v^)V0™||2 < ||A(/i™V0™)||2||V(/)™||ioe 

(7 19) 

< 4(||A//,„V(/.„J2 + j|VAi,„A(/),„||2 + ||/i™VA(/)„f )||V(/)™||i„. 
We treat the RHS of (|7.19p term by term. First we have 

(7.20) ||A^„V0„j2||V(/),„||ioo < i|A/i,„l|2||V0„Ji... 

By the Sobolev interpolation inequality ||/|||oo < C||/||L2j|/||^3 in 3D and (|3.2ip we have 

||A/i,„in|V0„||i,^ < C||AM™|P||V0„||i2||V0„||2,3 

(7.21) <C||A/i,„|P||V0„|p||A2(/,„,f 

<C||Am„.|P||A20„||2, 

where we have used the estimate of ||V(/)„i|p. Therefore, by plugging (|7.2ip into (|7.20p we have 

(7.22) ||A^„V0,„in|y^^^j|2^ <C||A/i,„in|A20,„||2. 
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For the second and third terms on the RHS of (j7.19p . using the same interpolation incquahty as above 
and Lemma 2.2 (X) and (XII) we have 

||VM™A<?i„|p||V0™|||oc, < ||VAim||i6||A<?i™||23||V0™||ioo 

< ||VAirn||ffi||A(/)„J||VA0„.||||V(/),„||||V(/)„jH3 

< ||AAt™|p||VA<?i„||||A20™|| 



(7.23) 



and 



(7.24) 



||MmVA0,„||2||V<?i™||i=e < ||Aim|Pl|VA<?i„||i^||V0™||||A20„|| 
<q|VA</)„,||||A30,„||||A2</,„J| 

< ^\\A^,„f + C||VA(/),„|n|A20,„|l2. 

Plugging (fr22|) - (fr24|) into (fTTOl we have 

(7.25) i|A(v„)V0„||2 < 1||A30„,||2 + C(||VA0„,||2 + ||Am„JP)||A20„||2. 
Plugging ((7Tf|) . ([7J5|) and t^TI^ into ((71^ we have 

(7.26) II A(v„, • V0„O||2 < i|| A30„||2 + C(||VA0„f + \\^„,\\l^ + || A^„f ) || A^^^f. 
Plugging (|7.26p into (|7.15p . after rearranging terms, we get 

(7.27) + W^'^rnW' < C{1 + l|VA0,„f + ||v,„ll2,, + |lAAt„j2)||A20,„f . 

Combining (fTTTI) and (fr27| we get 

|(||VA0„„ip + il A2</)„.iP) + (||VA20„||2 + ||A^0,„|p) 

< C(l + ||VA0„J2 ^ ||^^^j2^^ ^ jlAM,„|P)(|lVA0,„|l2 + |1A20,„||2). 

Step 5. Let T > be the life span of the solution. For any < s < t < T, applying Gronwall's 
inequality to (j7.28p over [s,t\ and using previous estimates we have 

(||VA0„||2 + ||A20„||2)(f) 
(7.29) - { / + + llvmilffi + \\Afi„,\\^){T)dT'j{\\'^^^mf + ||A20,„||2)(s) 

<C(||VA(/)„j2+||A20„j2)(s). 

Integrating (|7.29l) w.r.t. s over (0, t) we get 



(7.30) t{\\VA^m\\' + \\A^^mf)it)<C (||VA0™|p + ||A20„||2)(s)ds<C(r), Vte(0,r), 

Jo 

which implies that 

(7.31) (|jvA0„|p + ||A20„j2)(i)< Vie(0,f), 
where C(r) is an increasing function of T. 
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Step 6. For any < s < t < T, integrating (|7.28p over time and using (|7.3ip we have 

(|1VA0,„||2 ^ j|A20„j2)(t) + /* (||VA2</,,„i|2 + i|A-Vrn||')(T)dT 



(7.32) <(||VA,.|r + i|A^,.|r)(.)+C(T) / (l + i|VA0.„f + ||v-„||l,.+||A,„.iP)(.) ^^ 



<^ + + ||VA0„,||2 + ||v„,||l,. + ||A^„||2)(r)dr < ^^^^ 



s s 
In particular, we have 



(7.33) /* (||VA20„J2 + ||A30,„||2)(r)ds < ^^^^ 



s 

For any a > 0, muhiplying (|7.33p by s" then integrating the resuh w.r.t. s from to t we have 

t pt ft ^ci 



(7.34) / s" (||VA20„||2 + II A-Vm||') (T)dTds < C(T) / s^-irfs = C(r)-. 

Jo Js Jo 

Changing the order of integration we then have 

Ts" r(||VA20„||2 + ||A^0„,||2)(r)rfrrfs= f ( Ts^dsVuvA 
JO Js Ja \Jo / 



(7.35) 



</)„J^ + II A^0„ II ^)(r)dr 

r l!ji(||VA20„,||2 + ||A30„||2)(r)dr. 
Jo " + 1 



Plugging (|7.35p into (|7.34p we have 

(7.36) rr"+i(||VA20„||2 + ||A-V„||2)(r)dT<C(r)^ 







a 



for any a > and < t < T. This completes the proof of Theorem 1.4. 

8. Gevrey regularity in 2D & 3D 

In this section, we study the Gevrey (spatial) regularity of the solution, and show that the strong 
solutions obtained in Theorem 1.1-Theorem 1.2 belong to a Gevrey class of regularity defined below. 
Since the elements in the Gevrey class have high-mode coefficients that decay exponentially in wave 
number to zero, as a result of Gevrey regularity, one can show that the standard Galerkin method based 
on the basis (|2.ip ~ (|2.3p converges exponentially fast, see e.g. [TTl [TH HOI Ull HI] ■ Throughout this section 
we will assume that </>''(x) G iJ^(fi). It should be pointed out that in [TS] the authors studied Gevrey 
regularity for general nonlinear parabolic equations with analytic nonlinearity using the fact that If" is a 
Banach algebra for s > ^. But, Foias and Temam [16| get a better result for the Navier-Stokes equations 
when uq G H^, which is not an algebra in either 2D nor 3D, by taking advantage of the quadratic 
nonlinearity of the NSE. Here, one probably can follow those arguments and reduce the condition on the 
initial condition. However, this might be very involved since our nonlinearity contains (p^. So, we follow 
|15j instead to simplify the presentation and we do not claim that we have the optimal result regarding 
the assumption on the initial data. It is also important to note that the eigenfunctions p.ip - (j2.3p span a 
subspacc of the periodic functions on [0, 2L] x [0, 21] x [0, 2h] which is invariant under the flow generated 
by (jl.ip . Due to this symmetry, the boundary conditions do not explicitly enter the analysis and we can 
proceed as if (jl.ip were endowed with periodic boundary conditions. 
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For the domain fl under consideration and t > 0, wc define the Gevrey classes G(^^(f2) = L)(A^/^e*"^^'^), 
where A = {I — A). These spaces are Hilbert spaces with respect to the inner products 

(8.1) Ma^^aj =11 + bf re^*(^+l^-|^)^^^ 
with the corresponding norms 

(8.2) ||.||^./.(,,) = ( E \v,?{l + \3?Ye^'^'^\^\'''''\'\ 

where the Vj and Wj are the Fourier coefficients of v and w, respectively, with respect to the basis in 
(j2.ip - (|2.3p where the functions cos k9 and sin k9 arc represented in terms of the exponential function e**^^ 
and its conjugate e"**"'^, c.f. [15] . 

We have the following lemma (c.f. [15]). 

Lemma 8.1. Forp>2 and t > 0, G^/^{n) is a Banach algebra. That is 



for any u,v € G'f^^(r2). 



Concerning Gevrey class regularity, we refer the readers to [16] for the solutions of the Navier-Stokes 
equations and to [15] for solutions of general nonlinear parabolic equations and to [27] [29] for the Benard 
convection in a porous medium. 

Now we turn to the proof of Theorem 1.5. We first rewrite the equations p.23p in terms of A as 

+ P,„(Vm ■ V(/)„,.) ^ -Apim + Mm, 

\ Mm = Pm [(0m + 4>f - {<Prn + 4>)^] - 2(/'™ + Acj^m, 
Vjri = PmPa-(MmV(/),„). 

Applying A^e*^^^^ to (|8.3P i and taking the inner product of the resulting equation with A'^e*^^''^ (jjm 
we have 

[A e ——,Ae (pm]+[Ae (v,„ • V(/),„), A e (p,n) 
(8.4) \ dt J 

= - (A2e*^^^^AM™,A2e*^^^Vm) + (A2e*^''Vm,A2e*^^'Vm). 

After integration by parts and using (18.31) 9 we have 

(8-5) = - {A^e'^''\{<p„. + 4>f - 30™), A3e*^''V„0 - \\A^e'^''''M? 

+ {A^e'^^'^ad^rr. + 0)^), A^e'^^'Vm) - 2p2e*^^^Vmf , 
which implies that 

1 ^ II x2 tA^/^ / ||2 I II i3 tA^f^ I ||2 , oil /l2 tA^^^ , ||2 

(8.6) 2dt 

<ePV^''V™||2 + a(||Ae*^^''(v™ . V0™)||2 + \\A\'^''\{4>„. + 0)3||2 + p2gM^/^^^^|j2^ ^ 



where we have used Poincare's inequality to A^/^e*'^ </>,„. Choosing e = 1/2 we have 
d_ 
di 



•7) 4ll0m|lG2 + ll0m|lG3 + 4||0,„||^2 < C(||v„ • V0„J^i + || ((/),„ + 0)^11^2 + ll^mllcO' 

/Vf t t t t t t 
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For the second term on the RHS of (j8.7p . since is a Banach algebra, we have 

(8.8) ll('/'m + ^)'llG? <C|l('/''n + ^)llS? <C(ll'/'™llSf + 

Now we deal with the first term on the RHS of (|8.7p . Since Gl is a Banach algebra, we have 
|lv,„ • < ||v,„||gi||(/),„||L/2 

< C||/l,„||gi||(?:),„|lL/2||0m||G? 

(8.9) t ^ (^t t 

< C(||(0,„ +0)3||^j + + U,n\\l.)U,n\\%2 

< C{U,n\\%2 + 101^6^* + U,n\\l2)UrnfG2, 

where we have used Poincare's inequality in various places. 
Plugging dHTHl) and ([H^ into we have 

(8-10) ||I<^™IIg? + II-^^IIg? < CiU^\\%. + m'e'' + ||0™||^.)(||0™rG. + l) ^ P{U„,\\l^;t), 
which implies that 

(8.11) Um{;t)\\k+ f Un.{;s)fa.ds< f P{\\<l>,n{-, s)\\U ^)ds + Ulfci- 

Jo Jq ' 

Let us consider the first term on the LHS of (jSlTI) . Since \\(l>°ra\\Q2 = P^0°J1^2 < Co||(/)°||^4, by a 
continuity argument we know that there exists a time T** > such that 

(8.12) ||0„(.,t)||^. <2||<J|^. + l<2Co||0°||l,. + l, Vte [0,T„7]. 
Since P{p; s) is an increasing function for p > and in s, from (|8.1ip and (|8.12p we have 

(8.13) U„,{-,t)\\l,<tP{2CoUYm + l-,t)+GoU°\\%., Vi£ [0,Tr]. 
Let T** be some small value of t such that 

(8.14) T**P(2Co||0°||l,4 + 1;T**) < GoU^'Wl. + L 
Then it holds that 

(8.15) |10„(.,t)||^. <2Co|l0"|P^4 + l, yte[0,T**], Vm>0. 

Upon passing to the limit with rn — > oo we conclude that (|8.15|) holds for the limit function. Moreover, 
from (|7.28p we know that 4>{t,x) e H*{il) globally in 2D and locally in 3D if the initial data belong to 
H'^{Q). Therefore, by repeating the process we know that the solution possesses global Gevrey regularity 
in 2D and local Gevrey regularity in 3D. This completes the proof of Theorem 1.5. 

9. Long-time dynamics 

In this section, we investigate the long-time behavior of the solution to the original problem. First, 
we show that cj) converges exponentially fast to (f> in the norm as time goes to infinity, provided that 
(f> lies outside the spinodal region where F"{(j)) ~ icfp — 1 > (recall F(0) — \{(fP' — 1)^) and the initial 
perturbation \\(jp — 0||^2 is sufficiently small. The proof is in the spirit of [J. It should be pointed 
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out that the idea for proof has been used in [3] to study the long-time behavior of strong sohitions to a 
coupled Cahn-Hilliard-Navier-Stokes system: 

11 = F'{<i>) - A,/), 



on bounded domains in 2D and 3D. Similar results are obtained therein, i.e., converges exponentially 
fast to in the norm provided that lies outside the spinodal region where F"{-) > and 110" — 0||^2 
is sufficiently small. We briefly explain the idea as follows. Since (f> is assumed to be strictly greater than 
■\/3/3, we note that F"[(f)) ~ 3(f)^ — 1 > on a small neighborhood ~ [4> ~ S, (j) + 5], where S > 
is sufficiently small. The idea is to solve a modified problem with F replaced by an auxiliary function 
F^ whose second order derivative is non- negative on R and coincides with F" {■) on I^. Then, under 
the smallness assumption on the initial perturbations, it can be shown that the solution to the modified 
problem is indeed the solution to the original problem and converges exponentially to the constant state 
as time goes to infinity. The result holds true for both the 2D and 3D cases. 

On the other hand, we show that, in 2D, without the smallness assumption on the initial perturbation, 
4> still converges exponentially fast to 4> in the norm as time goes to infinity, provided that L, the 
length of the longest edge of the rectangle or the box, is strictly less than tt. This condition will trigger 
a chain reaction leading the energy estimate performed in Section 4 to a whole new scenario. The result 
is in strong contrast with the one above. Physically, the result indicates that making the thickness of 
the diffuse interface relatively large (compared to the dimensions of the domain) leads to constant phase 
states. 

In both situations, the key ingredient of the proofs is the estimate (|3.24p derived in Section 3. 

Remark 9.1. In this section, C and d will denote generic constants which are independent of the 
unknown functions and time. 

Step 1. We begin with some preparations. First, let us recall 



< 



Vf + (v • V)v + VP = ;^Av + 



V • V ^ 







(9.1) 



F{ci^) = -{ci^^-lf. 



Then it is obvious that 



(9.2) 



F' (</))= 03 - 0, F"((/)) = 302 - 1 and observe that a* = ^'(0) - A0. 



Second, we let 



(9.3) 




and assume that 



(9.4) 




Third, for fixed 



(9.5) 
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we set 

(9.6) h = [v-5,v + 5]. 
Then it is easy to see that 

(9.7) F"{(j)) = - 1 > on Ig. 
Fourth, one can easily construct a function F^, of class such that 

(9.8) F,(z.)=0, Fl{v)^Q, Fl'iz)^F"{z) on Is, ^^^(z) > if z i Is, 
and such that i^"' is bounded on R. By construction, it is easy to see that 

(9.9) F, > 0. 
By the convexity of F^ we see that 

(9.10) F,{iy)-F,iz)>F:,{z)iiy~z), 
which, together with (|9.8p . implies that 

(9.11) F^{z) < Fl{z){z^iy), V z G M. 
Moreover, there exist positive constants Fi and F2 such that 

(9.12) \Fiiz)\ < F, + F2\z\^ |F;(z)| <Fi+F2|zp, V z € R. 

Step 2. Next, we solve the original problem with F replaced by F^ in (IOJ2 and with the same initial 
data, and denote the solution by (pvj^^v and v^. In this situation, the estimate p. lip is 



(9.13) -(^ j^F,(0,)dx+-|jV0,||2 ) +||VAi.i|^ + l|v.r = 
Integrating (j9.13p in time we have 

(9.14) / F,{Mm^+\\\^MtW+ / {\\y^lu{s)r + \\Ms)f)ds< / FM'')d^ + ^\\^4Pf 



Using (|9lT]) and (|9?T2l) we have 

FAf)d^< [ F:(0"(x))(0"(x)-i.)dx 

<||f;(0°)||||0°-HI 
(9-15) < ||Fi+F2|0"|3||||0"-H| 

<C{l + \\ni^)U'-H 
<C{l + \\n'm)\\f-H- 

Plugging (|9.15p into (|9.14p we then have 



(9.16) 



1, 



F,(0.(t))rfx+-||V0,(t)||2 + (||VM.(s)|l2 + ||v,(s)||2)ds 
<C(1 + ||0O||?,OII0°-HI + ^1|V0"|P. 



Step 3. In this step, we exploit the estimate (|9.16p to get more information about the solution. Here 
we present two different versions of the proof. 
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Proof 1. Taking the L'^ inner product of with — v and noting that {fli^,(/)^ — v) = 0, where 
P-'^ ^ jn\ In t)dx, we have 

(9.17) {fi.-flu,(jju-iy) = \\V^,f+ [ FU(b,){(l),~i^)d^>\\V(j),M^+ [ F,{4>,)d^, 

Jn Jn 

where we have used (|9.1ip . 

Using Poincare's inequahty we have 

(9.18) (Ai. - -'y)<y.- fi.WU. - < c\\v^IMVM\ < c\\v^l,\\^ + ^W^^Ml^ 

which, together with (|9.17p . yields 

(9.19) IIVm.II' >C||V0,|p + C / F,i^,)dx. 



Taking the inner product of iii, with —Aip^, integrating by parts and using the relevant boundary 
condition we have 

(9.20) ifi,, = |lA0,f + / F;((/.,)|V0,pdx > |lA(/.,f , 
where we have used (j9.8p . By Poincare's inequality we have 

(9.21) (/i,, -Acj),) = {Vfi,,V(b,) < ||V/i,||||V0,|| < C\\Vfi,f + i||A0,||2. 
Combining (|9^ with (p:2T1) we get 

(9.22) llV^i^f >qiA0,f. 
Therefore, plugging (^1^ and ((O^ into ^^JE^ we have 

F,(,/),(O)dx+i||V0,(t)f+ 

(9.23) (c(||Va.,(s)||2 + ||A,^„(.s)f + llV0,(,s)f + ^ i^,(0,(s))dx) + |lv.(.)f )ds 

<C{l + UYm)U"-iy\\ + ^\\Wcl>Y, vt>o. 
Proof 2. Recalling the estimate p.24p we have 

(9.24) - + ||A0,J2 = _ ^ F"{cf>,)\W^.\'d^. 
Since F[,' > on M, we have 

(9.25) 1^110^ _^||2 + ||A0,||2<o, 
which gives 

(9.26) ^WMt) -'^f+ f Wu{s)fds < - HI'. 



Combining (|9.16p and (|9.26p we have 
(9.27) 



<c(i + UYh^)U" hi + ^l|v<^°|p + - i^r, V t > 0. 
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9.1. Long-time behavior in 2D for small initial perturbations. We first show the long time dy- 
namics in 2D. Since the solution exists for all time in 2D, when \\(fP — i^||^2 is sufficiently small, thanks 
to (IPs)) or (IOTI we have 

(9.28) / {\\V^l,{s)f + \\AMs)f)ds<ei, V t > 0. 

"'0 

By (p:^ we have 

(9.29) ||VA<^,||2<C(||VM.f + ||A0,f), 
which, together with (|9.28p . implies that 

(9.30) / \\WA(t)^{s)fds <e2, V t > 0. 

Jo 

By (|4.5p we have 

(9.31) ^l"^'^''"' + ^ll^''^-!!' ^ CWWAM^Acb.r + C||VA0,||2, 

where the constant C is independent of time. Applying Gronwall's inequality to (j9.3ip we have 

(9.32) !|A0,(i)i|2 < cxpjc^ ||VA0,(s)||2rfs|(^j|A/||2 + C^ \\VA(j>,fds^<e3, V i > 0, 

where we have used (|9.30p and the smallness of 110° — v\\'\j2- Therefore, by p.2ip we know that 

(9.33) \\4>.{t)-y\\H2 <£4, Vi>0, 
which, together with the Sobolev embedding ^ C°, implies that 

(9.34) \\(j)^{t) - u\\lo <e5, Vt>0. 

We remark that the constants Si ii — 1, • • • , 5) are independent of time and go to zero as \\(ff^ — ^'||^2 
tends to zero. 

We observe that, when \\(tP — i^|j^2 is sufficiently small, we can guarantee by (|9.5p that 0^ €E Is for any 
t>0 when we choose £5 < (5. Thus, 

(9.35) F'J{(p,) = F"{cj,,), 
which, together with (|9.8p . implies that 

(9.36) i^,(0,) = F(0,)-F(i/)-(0, -z/)F'(i/), e/i. 

From the definition of the solution we see that adding the above affinc function to F does not change the 
equations. Therefore, we conclude that (f>^ is a solution of the original problem. By uniqueness we know 
that — (j) for any t > when ||(/)'^ — v\\'\j2 is sufficiently small. As a by-product, we know that (f) E Is 
for any t > 0. 

We note that 4> satisfies the same estimate (|9.24p as 4>i, does: 

(9.37) - + II A0IP = - ^ F"(,/.)|V0pdx. 
Since F" > on and (fi G Is for any i > 0, we have 

(9.38) 1^110 _^||2 + II A0||2<o, Vi>0. 
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Poincare's inequality (c.f. (|3.2ip ) then gives 

(9.39) l|||0_^||2 + c|l</)-HP<O, yt>o, 
from which wc derive 

(9.40) U-i^lWt) <\\(^'> -u\\^e-^\ Vt>0. 

Since ||0 — i^||^2 is uniformly bounded by virtue of (|9.33p . by the interpolation inequality < 
C\\D^f WW f\\ for V/ • n|9o = we have 

(9.41) ||^-Hll,i(i)<Ce-^*- 

9.2. Global existence and long-time behavior in 3D for small initial perturbations. Now we 

turn to the long-time behavior of the solution in 3D. First, we need to show that the solution exists for 
all time when the initial perturbation is sufficiently small. We still work on the auxiliary problem with 
F replaced by F^. 

For this purpose, we re- visit the estimate (|5.14p . which reads 

(9.42) ^ll^'^'-ll' + W^^'^-W ^ CiWAc^4'' + C2|1A(/.,|12, 

for some constants Ci,C2 > independent of time. We also have the following estimate from (|9.23p or 
(im) : 

(9.43) / ||A0,(s)||2d5<C3||0°- 1/11^2, Vi>0, 

Jo 

for some constant C3 which is independent of time. Let 



(9.44) 110° - < e, where e<min|i 



4' 16(Ci(i)8 + C2)C3, 



and 



(9.45) r = sup{T>0: sup ||A(/)^(i)||2 < 1/2}. 

0<t<r 

Then it is clear by the continuity argument that T > 0. If T = -l-oo, then we have a global solution in 
3D. Suppose that T < +00. Let us focus our discussion over the interval [0, T). From (|9.42p we have 

-^||A(/).|p + ||A20,||2 < (Ci||A0,||i« + C2)||A0,|p 
(9-46) / 1 8 \ 

- (^1(2) +^2jl!A0.f. 

Integrating (|9.46l) w.r.t. t we have 

II A0.(t)|p < (c, Q\ C2) W^Ms)rds + WAcj^T 



(9.47) 











< 




[^'- 


FC2 




3 


1 


1 


< 










I6 + 


16 ~ 


4' 



iO- 1.11^2 + 110° 



This implies that 

(9.48) sup ||A0,(t)||2<i, 

0<t<T 4 
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which is a contradiction to the definition of T. Therefore, it must holds that T = +00. We conclude that 
if the initial data satisfies (|9.44p . then || A0y(t)|p < 1/2 for all t>Q. In particular, the solution exists 
globally in time in the 3D case. Moreover, we have from the above 

(9.49) \\l^Mt)f <(cA]-)\o^C:ie + e^Cie, V i > 0. 



From (|9.49p and Sobolev embedding we know that when e is sufficiently small, 110!^ — i^Hioo is sufficiently 
small for all time. Therefore, will be confined in the interval Is for all time. Since (|9.37p is valid for 
3D also, following the arguments in preceding section we know that the decay estimate still holds in 3D. 
In other words, we have 

(9.50) ||0(t) - z/||^i < Ce"'^*, Vi>0. 

9.3. Long-time behavior in 2D for large initial perturbations. In this subsection, we study a 
different aspect of the long-time dynamics of the solution in 2D. As pointed out in the beginning of this 
section, we show that, in 2D, without the smallness assumption on the initial perturbation, the solution 
4> still converges exponentially fast to 0, as time goes to infinity, provided that, as explained in the 
Introduction, the longest edge of the domain is shorter than tt. Moreover, the decay estimate will be 
extended to the norm of the solution. We begin with the estimate (|9.37|) : 

1 d 



In 



2 dt' 

(9.51) = - / (302 - l)|V0pdx 

n 



<||V0f <^||A0f, 



which gives 



(9.52) l|||<^_^||2+(i_g)||A0j|2<o. 
Since L < tt, integrating (j9.52p over time we have 



/ ||A0(s)|pds < C, Vt>0. 
Jq 



(9.53) 

Since (c.f. (1X131) ) 

(9.54) f \\\'^l{s)fds <C, Vt>0, 

Jo 

together with ((3?29)) . (|933)) and ((9?54l) we have 



(9.55) 

Recalling (|4.5p we have 



/ \\VA(f){s)fds <c, yt>o. 
Jo 



(9.56) ^I'l^'^ll' + ^ll^^'^ll' ^ C^l|VA0f j| A0f + C||VA0||2. 
Gronwall's inequality together with ()9.53p and (|9.55p then implies that 

(9.57) \\A(j){t)\\'^+ f \\^^(j>{s)\\'^ds<C, Vi>0. 
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Plugging the uniform estimate of ||A0jp into (|9.56p we have 

1 II A ,ll9 1 

where we have used (j3.29p . 

Next, we shall re- visit the estimate (|3.11l) which reads 



(9.58) oT^^n' + 7:\\^'n' < C||VA0|p < C,{\\S/fA\' + |lA0f ) 



(9.59) 
We note that since 

(9.60) F(0) = Fii^) + F\v){<i> + i.g(0)(0 - where g[<i>) = \{2v' + (0 + vf -2)>-\ 
and 

(9.61) / (</) - i^)dx = 0, 
it holds that 

(9.62) ^ / F(0)dx=i^ / g[(l)){<^-vfd^. 



dt Ja 2dt 
Plugging (p:^^ into (P^^)) we have 



(9.63) 



I Q g(</))(</. - ^)2dx + i|| V0|p) + II V/if + ||v|p = 0. 



Combining (|9.52p and (j9.63p and dropping the non- negative term ||v|p we have 

d /I 
di[2. 



(9-64) id I 5(0)(0-i^)'rfx+i||V</)||2 + l||0-HP)+||VM|P + (l-iV'r')llA0r <0. 



Let 



(9.65) E{t)^\ I g{^){cj,~vfdi,+h^<f>r+\u-vr>-\ I {cj>-vfd^+]-\\^<j,r+\u-v 

where we used (|9.60p . According to Poincare's inequality we have 

(9.66) E{t)>\(l-^)\\V<l^f + \u-vr. 
By absorbing the RHS of (^35)) into the LHS of (p:M|) we have 

Jt ( i_^^%2 ^W + ^l|A0|p) + CsdlV/if + ||A<^|p) + i|| A^0|p < 0. 
After integrating in time we have 

(t^^^^*^ + ^l|A0(t)f ) + £ (c5(||VMs)||^ + ||A0(s)|p) + i||A2<^(s)|p)ds 



(9.68) 



In particular, we have 

(9.69) ||A0(t)||2 < C, Vi>0, 
which implies that 

(9.70) ||(/>(t) - 1^11^2 < C, and <C, V t > 0. 
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Thanks to (|9.60p and the definition of the function g{(j)) we have 



1 

(9.71) 2 Jo 



5(0)(0-j/)^dx< 



2 



2 

2 



1 



110- HP 



<cU-'^\\, 

where we used (|9.70p in the last inequality. Therefore, by virtue of (|9.7ip and (j9.66p we see that 

(9.72) E{t)^U-uf + 

where ^ stands for the equivalence of quantities up to a multiplication by a constant. 
Let 

(9.73) K{t) ^ ^-^E{t) + Qii) ^ C^5(||VMf + || A0f ) + -!| A^-^lp. 

Then from p.2ip we know that there exists a constant Cg > such that 

(9.74) C^K[t) < Q{t), 
which, together with (|9.67p . yields 

(9.75) j^K{t) + CeKit)<0. 

Thus, K{t) decays exponentially to zero as time tends to infinity. By (j9.73p we then have 

(9.76) \\cf>^,^\\l.{t)<Ce-^'. 

Remark 9.2. The decay estimate obtained in this section can be extended to higher order norms of the 
solution in 2D, which can not be achieved by the method used in preceding sections. The penalty here is 
that the size of the domain can not be arbitrary. 

10. Conclusion 

In this paper, we studied the question of well-posedness and long-time qualitative behavior of a mixture 
model for solid tumor growth. Specifically, we have shown that: 1) when the initial data belong to H^{Vl), 
strong solutions (see Definition 2.1) to the IBVP are globally (locally resp.) well-posed in 2D (3D resp.); 
2) the spatial regularity of the solutions automatically increases by degree 2 within the lifespan of the 
solutions; 3) when the initial data belong to H^{Vl), the solutions are globally (locally resp.) spatially 
analytic in 2D (3D resp.); 4) in both 2D and 3D, the scalar volume fraction converges exponentially 
rapidly to its spatial average over the domain, as time goes to infinity, provided that the solution lies 
outside the spinodal region and the initial perturbation is sufficiently small; and 5) in 2D, if the length of 
the longest edge of the domain is strictly less than a constant multiple of the interface thickness, then the 
volume fraction still converges to its spatial average over the domain exponentially rapidly, as time goes 
to infinity, even if the amplitude of initial perturbation is large. The long-time behavior of the solutions 
suggest that the distinction between the tumor and the surrounding tissue in the microenvironment will 
blur and homogenize as time proceeds. On the other hand, it is unknown what happens if the solutions 
are in the spinodal region. This launches a new interesting problem for future pursue. In addition, we 
expect to extend the results to the case with non-divergence free velocity fields by adopting similar ideas 
to the ones presented in this paper. 
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